Pattern formation and glassy phase in the (f)^ theory with screened electrostatic 

repulsion 
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We study analytically the structural properties of a system with a short-range attraction and a 
competing long-range screened repulsion. This model contains the essential features of the effective 
interaction potential among charged colloids in polymeric solutions and provides novel insights on 
the equilibrium phase diagram of these systems. Within the self-consistent Hartree approximation 
and by using a replica approach, we show that varying the parameters of the repulsive potential 
and the temperature yields a phase coexistence, a lamellar and a glassy phase. Our results strongly 
suggest that the cluster phase observed in charged colloids might be the signature of an underlying 
equilibrium lamellar phase, hidden on experimental time scales. 
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Introduction — In the last few years there has been 
much interest in the role of the inter-particle potential 
on controlling the structure and the dynamics of colloidal 
suspensions [llSSSIal due to the potential apphcation 
of these systems for designing new materials with a wide 
range of viscoelastic properties. In charged colloidal sys- 
tems the effective interaction can be described in terms of 
a short-range attraction and a long-range screened elec- 
trostatic repulsion (well approximated by the DLVO po- 
tential 0). The competition between attractive and re- 
pulsive interactions on different length scales stabilizes 
the formation of aggregates of an optimal size (cluster 
phase) characterized by a peak of the structure factor 
around a typical wave vector, km, that has been ex- 
perimentally P, and numerically H, Q |^ observed. 
By an appropriate tuning of the control parameters, the 
system progressively evolves toward an arrested gel-like 
disordered state {colloidal gelation). Although intensely 
studied both experimentally and numerically, a theoret- 
ical understanding of the mechanisms underlying these 
phenomena is still lacking and many gaps remain in our 
knowledge of the equilibrium phase diagram of these sys- 
tems. 

In this letter we study analytically a i/)^ model with 
competition between a short-range attraction, described 
by the Ginzburg-Landau Hamiltonian, and a long-range 
screened repulsion, described by a Yukawa potential. Al- 
beit schematically, this model contains the essential fea- 
tures of the effective interaction among charged colloids 
and sheds new light on the structural properties of these 
systems. Depending on the control parameters, there is a 
region of the phase diagram where usual phase separation 
takes place. Conversely, as the screening length and/or 
the strength of the repulsion exceeds a threshold value, 
phase separation is prevented. In this case, at moder- 
ately high temperature the competition between attrac- 
tion and repulsion has the effect to produce modulated 



structures, which in the terminology of particle systems, 
correspond to the cluster phase. These modulated struc- 
tures are the precursors of a first order transition towards 
an equilibrium lamellar phase found at lower tempera- 
tures. By using a replica approach for systems without 
quenched disorder 0, and employirig the Self-Consistent 
Screening Approximation (SCSA) j^], we also show the 
presence of a glass transition line in the low tempera- 
ture region, once the first order transition to the lamellar 
phase is avoided. Note that the mechanism for the glass 
transition in this case is not due to the presence of hard- 
core type of potential. Instead, it is due to the formation 
of the modulated structures which order up to the size of 
correlation length. The geometric frustration, resulting 
in arranging such modulated structures in a disordered 
fashion, leads to a complex free energy landscape and, 
consequently, to a dynamical slowing down. 

Our results suggest that the cluster phase observed in 
colloidal suspensions should be followed, upon decreas- 
ing the temperature (or increasing the volume fraction), 
by an equilibrium periodic phase (a tubular or a lamellar 
phase, depending on the volume fraction). If, instead, 
this ordered phase is avoided, a structural arrest, corre- 
sponding to the gel phase observed in the experiments 
and in numerical simulations, should eventually occur. 
The existence of ordered phases in colloidal suspension 
and the fact that the transition to the gel phase could 
occur in a metastable liquid, are novel predictions, which 
have never been considered before. Recently, motivated 
by our work, de Candia et al. have unambiguously 
shown the presence of such ordered phases by using MD 
simulations on an atomistic model system of charged col- 
loids, interacting via the DLVO potential. 

Model and phase diagram — We consider the standard 
three dimensional <p'^ field-theory with the addition of a 
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repulsive long-range Yukawa potential: 
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where ^(x) is the scalar order parameter field. The model 
has been also studied in |0| in the context of microemul- 
sion. The parameters W and A are, respectively, the 
strength and the range of the repulsive potential. For 
W = we obtain the canonical short-range ferromag- 
net. Interestingly, for A — > oo we recover the case of 
the Coulombic interaction [H El El El El El El- 
This model has been used to describe the phenomenol- 
ogy of a wide variety of systems, where competing inter- 
actions on different length scales stabilize pattern forma- 
tions and the creation of spatial inhomogeneities (for a 
review see jl^ ) • These systems include magnetic systems 
and dipolar fluids characterized by long-range Coulom- 
bic interactions Ell j mixtures of block copolymers [H , 
water-oil-surfactant mixtures and doped Mott insu- 
lator, including the high Tc superconductors m. As a 
consequence, our model allows to describe and to inter- 
pret in an unified fashion the phenomenology of a wide 
variety of different systems. 

Here we present only the main results, and we refer to 
a paper in preparation for more details m . In Fig. 
the phase diagram of the model as function of the tem- 
perature, T, and the strength of the repulsion, W , is pre- 
sented, for a fixed value of the screening length, A = 2. 
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FIG. 1: (color online) Temperature (r)-repulsion (AtiW) 
phase diagram of the model for A = 2 (and ro = —1), showing 
the relative positions of the paramagnetic (P), ferromagnetic 
(F), lamellar (L) and glassy (G) phases. The value of g is such 
that TciW = 0) = 1. The continuous (black) and the dashed 
(blue) curves, found within the Hartree approximation, corre- 
sponds respectively to the second-order phase transition from 
P to F, TciWyX), and to the first-order transition from P 
to L, Tl{W,X). The dotted (green) and the dashed-dotted 
(red) curves, found within the SCSA, identify, the dynamical, 
TdiW, A), and the ideal, Tk(W, A), transition temperatures to 
G. The thin dotted (black) line (^ — Im) marks a crossover 
temperature below which the system establishes modulated 
structures, schematically sketched in the figure. 



Ferromagnetic and lamellar phases — We first solve the 
model within the self-consistent Hartree approximation, 
which consists in replacing the term g(t>^ /A of Eq. |^ with 
3(0^)0^/2 24]. This substitution allows to compute the 
correlation function G(k) = (^ki/'-k) — (0k) (0-k)- In the 
paramagnetic phase, (0k) = 0, one obtains: 

rG"i(k)==r + fc2-H47rIF/(A-2 + fc2)^ (2) 

where the renormalized mass term, r, is defined as: r = 
ro + 3(7(02). Since (0^) = ^^^^^^ ^ G(k), from Eq. @ 
the following self-consistent equation for r is derived: 

(A is an ultraviolet cutoff). For convenience let us define: 

AttWc = A~l (4) 

For W < Wc and a fixed value of A, we find a line of 
ordinary second order critical points, Tc{W,X) (continu- 
ous curve in Fig.^, separating a high-temperature para- 
magnetic phase from a low-temperature ferromagnetic 
one. This transition is characterized by the divergence 



of the susceptibility, i.e. G^^(k = 0) = 0, with the usual 
Hartree critical exponents (e.g., 1^=1 and 7 = 2 in three 
dimensions). Thus, for W < Wc usual phase separation 
occurs: the effect of the repulsive interaction is to de- 
crease the value of Tc from the critical temperature of the 
standard Ginzburg-Landau model (for IF = 0) to zero 
temperature (for W ^ Wc). Conversely, above Wc there 
is no phase separation. This result is quite important 
for designing new materials as well as in the experimen- 
tal and numerical study of colloidal systems, where it is 
crucial to distinguish the slowing down due to colloidal 
gelation from that due to kinetic of phase separation. 
Interestingly, the threshold value Wc, Eq. Q), coincides 
with that estimated for an atomistic model system of 
charged colloids interacting via the DLVO potential |5|. 

In the limit, A 00, we recover the already known 
results for the Coulombic case: phase separation occurs 
only for IF = El El El El El El E3 • 

For W > Wc, the model instead exhibits a first order 
transition fine Tl{W,X) (dashed curve of Fig. sepa- 
rating the paramagnetic phase from a lamellar phase |l7l 
m, characterized by a spatially modulated order: 

(0k) = A(<5(k-k„0 + ^(k + k„0), (5) 
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with amplitude A and wave vector km, given by 



(6) 



For W > Wc the spinodal line of the "supercooled" para- 
magnetic phase is located at T = 0, where the suscepti- 
bility diverges {TG~^(km) = 0). As a result, the first- 
order transition to the lamellar phase can be kinetically 
avoided and l ong -time g lassy relaxations can be, instead, 
observed [Illillllllilillilii. 

Glass transition — In order to analyze the glass transi- 
tion in our model, we employ a replica approach formu- 
lated to deal with systems without quenched disorder ■ 

The equilibrium free energy, F — —TlnZ, is relevant 
only if the system is able to explore ergodically the entire 
phase space. This is not the case, of course, in the glassy 
phase, where the system is frozen in metastable states. 
In order to scan the locally stable field configurations, we 
introduce an appropriate symmetry breaking field V'(x), 
and compute the following partition function j^]: 

m ^Jvcb exp (^-pn[(f>] - ^/d'x [V'(x) - 0(x) 

where m > denotes the strength of the coupling. The 
free energy /[ip] = — TlnZ['i/;] will be low if equals 
to configurations which locally minimize 7i[0]. Thus, in 
order to scan all metastable states we have to sample all 
configurations of the field -0, weighted with exp (—/?/[■(/;]): 



F= lim 

11^0+ 



(8) 



_F is a weighted average of the free energy in the vari- 
ous metastable states; if there are only a few local min- 
ima, the limit behaves perturbatively and F equals the 
true free energy F. However, in case of the emergence of 
an exponentially large number of metastable states with 
large barriers between them, a nontrivial contribution 
arises from the above integral even in the limit u — s- 0^ 
and F differs from F. This allows to identify the com- 
plexity, S, via the relation F = F — TE 0. In order to 
get an explicit expression for S we introduce replicas: 



F(m) 



T 

lim — In 

u-*o+ m 



ViJ; {Z[il; 



from which, we obtain F = dmF{m) / dm\m=i and 



S = T 



(j)F{m) j dm\ 



(9) 



(10) 



Integrating Eq. @ over we get an action which is for- 
mally equivalent to the (replicated) action of a system in 
a quenched random field. We can thus use the SCSA 
a technique developed to deal with such systems, which 
allows to determine the correlators in the replica space. 
The SCSA amounts to introduce a A^-component ver- 
sion of the model and summing self-consistently all the 



diagrams of order 1/A^ [H, [3. Since the attractive 
coupling between replicas is symmetric with respect to 
the replica index, one can assume the following struc- 
ture of the correlators in the rephca space: Gob(k) = 
[G'(k) — i^(k)] ^ab + -F(k), i.e., with diagonal elements 
G(k) and off-diagonal elements -F(k). For systems with 
quenched disorder, this ansatz turns out to be equivalent 
to the one-step replica symmetry breaking. While the 
diagonal correlator can be interpreted as the usual one- 
time equilibrium correlation function, rG(k) = (0k0-k), 
the off-diagonal term can be interpreted as measuring the 
long-time correlations: TF(k) = limt^oo((/)k(i)(/'-k(0)). 
Hence, F(k) vanishes in the paramagnetic phase while is 
finite in the glassy one. 

The system undergoes a glass transition in the low tem- 
perature region for W > Wc, with exactly the same na- 
ture of that found in mean-field models for glass- formers. 
Lowering the temperature we first find a purely dynami- 
cal transition at temperature (dotted curve in Fig.^). 
Here, the complexity, Eq. (|10|) . jumps discontinuously 
from zero to a finite value, signaling the emergence of 
an exponentially large number of metastable state. The 
complexity decreases as the temperature is decreased and 
vanishes at Tk (dashed-dotted curve in figure) where the 
thermodynamical transition takes place. 
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FIG. 2: Main frame: Momentum dependence of the corre- 
lation function, G(k), for 4ttW = 0.2 and A = 2 at T = Ta, 
showing that it is peaked around the typical modulation wave 
vector km with broadening given by the inverse of the 
correlation length. Inset: Momentum dependence of the non 
ergodicity parameter /k for the same values oi W, X and T. 

Structural properties — The correlation function, G(k), 
is plotted in Fig.|21at the dynamical transition tempera- 
ture T(i, showing a maximum at k^, defined in Eq. 10, 
with width given by the inverse of the correlation 
length. The correlation function in the real space reads: 
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G(|x|) ~ e~l''l/^ sin(fcm|x|)/|x|. This expression implies 
that, ahhough no periodic order occurs ((0k™) = 0), a 
lamehar structure of wave length l„i = over a fi- 

nite range f is formed (as sketched in Fig. [T]). In the 
glassy phase, T < T^, one has that ^ > 21^ llg: thus 
these modulated structures form over a length larger than 
their modulation length and become frozen. The glass 
transition arises from the fact that there are many possi- 
ble configurations to arrange such modulated structures 
in a disordered fashion, leading to a great number of 
metastable states. The presence of this characteristic 
wave length dominates also the dynamics, as indicated 
by the momentum dependence of the non ergodicity pa- 

rameter, /k = Mmt^^ t^^lm-'ti^ = iS}' P^°"^^ 
the inset of Fig. |21 at T^. The presence of a maximum 
at signals the fact that structural arrest is more pro- 
nounced over length scales of order . At higher temper- 
atures, T > Td, the non ergodicity parameter, /k, van- 
ishes and the glassy phase disappears; correspondingly, 
G'(k) broadens and the height of the peak decreases; 
hence, ^ decreases until the modulated structures fade 
continuously, approximately at a crossover temperature 
where ^ Im (dotted curve in Fig. P). 

These results are intimately related to the phenomenol- 
ogy observed in colloidal systems, where the competi- 
tion between attraction and repulsion leads to the for- 
mation of aphase of stable clusters at low tempera- 
tures P, U 0j J which is the analog of the modu- 
lated structures here found. Our results suggest that the 
transition to the disordered gel phase, numerically and 
experimentally observed in colloidal suspensions, occurs, 
in fact, in a metastable liquid, due to the presence of an 
underlying equilibrium lamellar phase (which might be 
more easily detected by increasing the screening length 
A). This novel prediction has been confirmed by recent 
MD simulations of a model systems of charged colloids 
in 3(i j^; a clear indication of the presence of periodic 
phases was also numerically found in 2d |29i |. 

Conclusions — We have derived analytically the com- 
plete phase diagram of a model with competition between 
short-range attraction and long-range screened repulsion, 
which contains the essential features of the interaction 
potential of charged colloids. To our knowledge, this is 
the first theoretical investigations on these systems. Our 
predictions have been confirmed by recent numerical sim- 
ulations and may be also experimentally checked. 
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